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SUMMARY 


A procedure is described for obtaining minimum-mass designs of insulated 
composite structural panels. The panels are characterized as consisting of a 
structural layer and an insulation layer. The panels are loaded by a general 
set of inplane forces applied to the structural layer and a time-dependent tem- 
perature applied to the outer surface of the insulation layer. Temperature and 
stress histories in the panel are given by closed-form solutions , and optimiza- 
tion of the insulation and structural thicknesses is performed by nonlinear 
mathematical programming techniques. Design constraints are enforced at a 
finite number of discrete times over the time period of interest, and satisfac- 
tory results are obtained with a small number of times. 

The computerized procedure is intended for preliminary design calculations 
to evaluate materials for specified applications and to perform parameter stud- 
ies. A set of design calculations is presented to evaluate the efficiency of 
the following eight structural materials under combined heating and mechanical 
loads: graphite/polyimide (Gr/PI) , graphite/epoxy (Gr/E) , boron/aluminum 

(B/Al) , titanium (Ti) , aluminum (Al) , Rene^ 41 , carbon/carbon, and Lockalloy. 

A study is also performed to assess the effect on design mass of the intensity 
and duration of heating for the eight materials. Examinations of final designs 
indicate that, for sufficiently high mechanical loads, an optimum structure has 
a temperature response well below the recommended allowable temperature for the 
material. A comparison of these designs with those assumed to operate at the 
allowable temperatures indicates that significant mass savings may be attainable 
by the lower temperature operation. 


INTRODUCTION 

A current technical challenge in the design of high-speed flight vehicles 
is the development of optimum structural design concepts for combined heating 
and mechanical loading applications. A concept of interest for applications to 
atmospheric reentry vehicles is the insulated configuration in which the struc- 
ture is protected from heating by insulation placed between the airstream and 
the structure. The potential benefits of 'structural concepts composed of fiber- 
reinforced composite materials are currently being evaluated (ref. 1). At an 
early stage in the design of such thermal-structural concepts, it is beneficial 
to conduct selected design studies to (1) evaluate candidate structural and 
insulation materials to narrow the choice of materials to a small number of con- 
tenders and (2) obtain a preliminary estimate of the optimum combination of 
structural and insulation material. An emerging philosophy for conducting such 
studies is integrated thermal- structural analysis and design (refs. 2 to 4) in 
which all important interactions between the thermal and structural responses 
of the structure are properly accounted for at every stage in the design pro- 
cess. Presented herein is a design technique to define the best combination 
of insulation and structural thicknesses which minimizes the total mass under 
a specified set of inplane mechanical loads and a transient heating pulse. 



Previous investigators have addressed the minimum-mass design of insulated 
structures (refs. 5 to 8) , but consideration has been limited to metallic struc- 
tures, uniaxial loading, and a restricted class of heating pulses. These sim- 
plifications were needed since the classical minimization techniques required 
extensive manipulations of algebraic equations. 

In the present work the structure considered is either a metal layer or a 
composite laminate under a general set of inplane mechanical loads and heating 
pulse. This generalization is facilitated by the use of nonlinear mathematical 
programming techniques to generate the designs systematically. An important 
aspect of the present approach is the treatment of the requirement that time- 
dependent temperatures and stresses be less than specified allowable values for 
all appropriate values of time. This aspect has received relatively little 
attention. In reference 9, which deals with the optimization of an ablating 
heat shield, constraints on time- dependent temperatures and stresses are 
replaced by time-integrated averages of the original constraints. The inte- 
grals are carefully formulated to avoid constraint violations, but, as pointed 
out in reference 9, the integral representation of the constraints leads to a 
smoothing effect which causes temperatures and stresses to be somewhat insensi- 
tive to changes in the structural size parameters. The approach used in the 
present work is a simple time slicing in which the constraints are enforced at 
a sufficiently large number of discrete times to ensure satisfaction of all con- 
straints throughout the time period of interest. (Since the completion of this 
work, progress has been reported in ref. 10 on methods to track one or more most 
critical times and enforce constraints only at those times.) Analytical solu- 
tions are obtained for the temperature history and stresses in the structure. 

The mass minimizations are carried out by use of the general-purpose optimizer 
AESOP (ref. 11). 

The analytical design procedure is applied to two studies: (1) A compara- 

tive structural efficiency study of several structural materials suitable for 
high- temperature structural applications - graphite/polyimide, graphite/ epoxy, 
boron/aluminum, titanium, Rene’ 41 , aluminum, carbon/carbon, and Lockalloy; 
and (2) a study of the effect of intensity and duration of heating on the 
required mass for these materials. Effects of certain simplifying assumptions 
made in the analysis are also assessed for selected final designs. The purposes 
of this paper are to describe the analytical design procedure and to present 
the results of the studies carried out with the procedure. 


SYMBOLS 

The values are given in both SI and U.S. Customary Units. The calcula- 
tions were made in U.S. Customary Units. 

c heat capacity, J/kg-°C (Btu/lbm-°F) 

c-|,C 2 heat capacity of insulation and structural material, respectively, 
J/kg-°C (Btu/lbm-°F) 

E-| ,E 2 Young's modulus in fiber and transverse directions, respectively. Pa 
(lbf/in 2 ) 
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F 


Tsai-Wu failure criterion 


F-| = 1/X T + 1/X C 

F 2 = 1/Y T + VY C 

Fll = -l/X>nX c 

f 22 = -1/Yt y C 

I 

F 66 = vs 2 

G shear modulus. Pa (lbf/in 2 ) 

g constraint 

k thermal conductivity, W/m-°C (Btu/in-sec-°F) 

k-j ,k 2 thermal conductivity of insulation and structural material, respec- 

tively, W/m-°C (Btu/in-sec-°F) 

N LAYER number of layers in laminate 

NTIMES number of discrete times 

N x' N y' N xy stress resultants, kN/m (lbf/in.) 

N tx' N ty thermal forces, kN/m (lbf/in.) 

RT room temperature 

S shear strength. Pa (lbf/in 2 ) 

T temperature, K (°F) 

T a allowable temperature, K (°F) 

T e g outer surface temperature, K (°F) 

* T eqo' T peak' T eqf '1 parameters of temperature pulse (fig. 2) , K (°F) 

T peak' T f J 

« T* optimum operating temperature, K (°F) (fig. 3) 

t()/t 45 ,t 90 thickness of 0°, ±45°, and 90° plies, respectively, cm (in.) 
t] , t 2 insulation and structure thicknesses, respectively, cm (in.) 

V Von Mises 1 stress criterion. Pa (lbf/in 2 ) 

W mass, kg (lbm) 
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X T ,Xc tensile and compressive strengths in fiber direction, respectively. 

Pa (lbf/in 2 ) 

Y T ,Y C tensile and compressive strengths in transverse direction, respec- 
tively, Pa (lbf/in 2 ) 

x,y laminate coordinate directions 

z-j,Z 2 coordinates through depth of insulation and structure, respectively, 
cm (in.) 

a lr a 2 coefficients of thermal expansion in fiber and transverse directions, 
respectively, °C -1 (°F~^ ) 

3 edge-fixity coefficient 

9 ply angle, deg 

V-j ,\>2 Poisson's ratio for fiber and transverse directions, respectively 

p mass density, kg/m 3 (lbm/in 3 ) 

P-|,p 2 mass density of insulation and structural materials, respectively, 
kg/m 3 (lbm/in 3 ) 

a a allowable stress. Pa (lbf/in 2 ) 

0 lr a 2' a 12 stress components in laminar coordinate system. Pa (lbf/in 2 ) 

t time, sec 

T period of time that constraints are imposed? monitoring time, sec 


CONFIGURATION AND LOADS 

The configuration, as shown in figure 1 , consists of a layer of insulation 
attached to a structural layer. This configuration is assumed to be a localized 
region of a larger structure and is sufficiently small so that neither the 
mechanical loads nor the temperature vary significantly over the region. Also 
the region is assumed to be sufficiently far removed from restrained boundaries 
so that edge effects are negligible. The structural layer is subjected to a 
general set of inplane mechanical loads given by the stress resultants N x , Ny, 
and Nxy, which are assumed to be time independent. The structural layer is 
either a filamentary composite laminate or a metal layer. No mechanical loads 
are applied to the insulation layer. The outer surface of the insulation is 
subjected to time-dependent heating, which is characterized by a temperature 
T e g(T) . in the present work, the particular form of T e q(T) is the curve shown 
in figure 2. The shape of the curve (generally a third-order function of time 
for T = Tf and constant thereafter) is determined by the quantities T e q 0 , 

T peak' T peak' T f' anc ^ T eqf* 
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ANALYTICAL DESIGN PROCEDURE 


Thermal and Structural Analysis 

Analysis of the temperature history is based on one-dimensional heat trans- 
fer through the depth of the insulation-structure configuration. The outer sur- 
face of the insulation has a prescribed temperature, and the back face of the 
structural layer is assumed to be an adiabatic surface. The temperature history 
throughout the insulation and structure is obtained by an analytical solution 
described in appendix A. The assumption that the back face of the structural 
layer is an adiabatic surface is known to be conservative. Quantification of 
the effects of the adiabatic assumption is described in appendix B f where this 
assumption is relaxed by permitting radiation from the back face of the struc- 
tural layer, and temperatures for radiative and adiabatic back-face conditions 
are compared. 

The structural analysis is based on calculating stresses and strains at a 
point in the structure resulting from a set of applied forces N x , Ny, and Nxy 
and the temperature at that point. This approach neglects edge-boundary effects, 
especially those due to thermal stresses arising from restraint of thermal 
expansion; however, thermal stresses due to interlayer thermal property mismatch 
are included. The applied forces are assumed to be constant (i.e., time inde- 
pendent) . The structural analysis based on these assumptions is described in 
appendix A. Effects of neglecting restrained thermal expansion due to edge 
restraint and effects of time-dependent loading are evaluated in appendix B. 

In the analysis, linear representations of the temperature dependence of the 
mechanical properties are used. In principle, any known variation of mechanical 
property with temperature may be input. 


Design Objective and Constraints 

The objective of the present procedure is to determine the insulation 
thickness t-j and structural thickness t 2 for the configuration in figure 1 
such that the total mass is a minimum and the structural temperature and 
stresses do not exceed prescribed allowable values throughout a span of time 
t. The mass is given by 


W = Pi t-j + p 2 t 2 (1) 

where for a composite laminate with 0°, ±45°, and 90° plies 

t2 = + fc 45 + t 90 (2) 

and for a metal structure t 2 is the metal thickness. The temperature con- 
straint is 
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(0 5s T = T) 


(3) 


T(t) S T a 


where T is the temperature of the structural layer and T a is the allowable 
temperature of the structural layer . The value assigned to T a is somewhat 
arbitrary and represents the maximum service temperature of the material for 
structural applications. In the case of filamentary composites, values of T a 
are generally dictated by the integrity of the resin material. In this paper, 
values of T a are based on current usage of the various materials in high- 
temperature applications and recommended temperature ranges from the appropriate 
literature (refs. 12 to 15). 

The stress constraint for each layer of a filamentary composite laminate 
is expressed by a modified form of the Tsai-Wu failure criterion (ref. 16) 

2 2 2 

F = F-| a-| + F 2 O 2 + Fn ai + F 22 G 2 + F 66 a 12 = 1 (0 < x £ x) (4) 

The coefficients F-j , F 2 r Fn , F 22 » and F 55 in equation (4) are functions 

of the allowable stresses for the composite system and are all temperature 
dependent. The stress constraint for a metal layer is given by 


F = v/a a S 1 


(0 S T £ T) (5) 


where the Von Mises' stress criterion is 


V = 


/ 2 2 2 \V 2 

(0, +02-0^0 2 + 3 ar 12 j 


(6) 


Before describing the analytical design procedure, it is instructive to describe 
the trade-off involved in designing an insulated structure. As shown in fig- 
ure 3, the required mass of the load-carrying structure increases with operat- 
ing temperature of the structure, primarily because of degradation of mechanical 
properties. The required mass of the insulation decreases with increasing oper- 
ating temperature level. The total mass therefore attains a minimum for some 
temperature T*. This optimum temperature may be significantly less than the 
allowable temperature indicated by T a . It is shown subsequently that the dif- 
ference between T* and T a is principally a function of the mechanical load. 


Design Formulation 

The design problem is formulated and solved as a nonlinear mathematical 
programming problem. The objective function is given by equation (1). The 
design variables consist of t-j and either the set tg^ fc 45 f and tgg or 
t 2 • The constraints are expressed by equation (3) and either equation (4) 
or (5) . For convenience and in preparation for the discretization of the time 
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dependence of the constraints, equations (3) to (5) are written as nondimen- 
sional constraint functions 


91 (T) = T(T)/T a -1=0 (0 S T S T) (7) 


gj +1 (T) = Fj(T) - 1 (0 S T S T; j = 1 , 2, . . . NLAYER) (8) 

where F is given by equation (4) or (5) . In addition to the behavior con- 
straints in equations (7) and (8), minimum gage constraints are applied to the 
design variables. 


Treatment of Constraints on Time-Dependent Quantities 

In developing the present analytical design procedure, special attention 
was focused on the satisfaction of constraints on time-dependent temperatures 
and stresses (eqs. (7) and (8)). Some guidance was available from reference 9, 
which treated the optimum design of an ablating heat shield. In reference 9 
the constraints on time-dependent quantities were replaced by time-integrated 
averages of the original constraints. The representations were carefully for- 
mulated to exclude violations at any point in time but, as pointed out in ref- 
erence 9, had the disadvantage that temperatures and stresses were somewhat 
insensitive to structural size parameters because of the smoothing effect of 
the integral representation of the constraints. Early in the present work, 
alternative approaches were considered, including other integral representa- 
tions (ref. 17) and ideas for tracking one or more most critical times and 
enforcing constraints only at those times. The approach decided upon was a 
simple time slicing, whereby constraints were enforced at a sufficiently large 
number of discrete times in the time period of interest. A schematic represen- 
tation of the discrete-times approach for temperature constraints is shown in 
figure 4. 

In the discrete-times approach, the constraints of equations (7) and (8) 
are replaced by 


gj[ = T(Ti)/T a -ISO (i = 1 , NTIMES) (9) 

9NTIMES+i*j = F j( T i) “ 1 = 0 

(i = 1, NTIMES; j = 1 , 2, . . . NLAYER) (10) 

Thus the total number of constraints is 


(1 + NLAYER) (NTIMES) 
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Because the number of constraints is proportional to the number of discrete 
times NTIMES, it is of interest to determine how large NTIMES must be for 
typical calculations. The answer is problem dependent , but studies carried out 
and documented in reference 4 and in a subsequent section of this paper indicate 
that, for the types of insulated structures under the types of heating histories 
considered herein, only a few (three or four) discrete times are needed. 


Sizing Technique 

The computer program used to perform the optimizations is a general-purpose 
optimizer denoted AESOP (ref. 11) , which accounts for constraints using an exte- 
rior penalty function approach (ref. 17). The AESOP code contains several opti- 
mization algorithms that can be selected in various combinations by the user. 

The approach used here was a combination of adaptive search and pattern search. 
These algorithms are described in reference 11. 

APPLICATIONS AND DISCUSSION 
Overview 

In this section the design procedure is used to carry out design studies. 
First a preliminary study is performed to assess the minimum number of discrete 
times needed to obtain minimum-mass designs. Next the design procedure is 
applied to an efficiency study for eight structural materials over a wide range 
of loads and two heating conditions representative of high-speed flight of 
winged reentry vehicles (refs. 18 and 19). The materials in the study include 
graph ite/epoxy, graphite/polyimide , boron/aluminum, aluminum, titanium, 

Rene 41, carbon/carbon, and Lockalloy. A study is then carried out to assess 
the effect on minimum mass of heating duration and intensity for the eight 
materials. Finally results are presented to assess the effects of certain sim- 
plifying assumptions made in the analysis. 


Materials Used In Studies 

The eight structural materials selected for the studies are listed in 
table I along with averaged values used for thermal conductivity, specific 
heat capacity, and density, as well as the allowable temperatures used in the 
studies. Properties are also given for the insulation which is LI-900 reus- 
able surface insulation. Mechanical properties for the structural materials 
were obtained from references 4 and 12 to 15 and are listed in table II. Min- 
imum gages considered are given in table III. 


Discrete-Times Convergence Study 

The purpose of this section is to show approximately how many discrete 
times are needed for a typical design problem. 
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The study is carried out for a Gr/E panel with an outer surface temperature 
history (fig* 2) characterized by Tp ea fc = 1089 K (1500° F) f Tpeak = 700 see, 

Tf = 1500 sec, and T e qf = 297 K (75* F) and loads of N x = Ny = -700 kN/m 
(-4000 Ibf/in.) and Nxy = 525 kN/m (3000 lbf/in.) • The monitoring time f 
was 6000 sec. Designs were obtained for 3 to 1 5 time slices. Final mass is 
plotted against the number of slices in figure 5. The plot indicates that only 
three or four discrete times are necessary to give satisfactory results in the 
sense that increasing the number of discrete times does not improve the result. 
Figure 6(a) shows the temperature and stress response (in the 45° ply) in the 
final design. The response of temperature is most important to this design 
(i.e., only the temperature constraint is critical for T > 0) . Observe that 
the temperature response is a smooth curve with a single local maximum. Thus 
only a small number of equally spaced time slices are needed to define properly 
the critical temperature. If a more severe or complicated heating pulse is 
applied, the temperature response would still be smooth and gradual. This 
assertion is based on previous analyses (refs. 4 and 20) and occurs because 
the insulation has a smoothing (dampinglike) effect on the thermal response of 
the structure. 

The stress-response curve is somewhat more complicated than the temperature 
response but again only has a single local maximum. Since stress is not criti- 
cal for t > 0, it is not clear from this example what effect the stress 
response has on the number of discrete times required. To investigate this 
effect, consider the Gr/E panel with higher loads of N x = Ny = -1.05 MN/m 
(-6000 lbf/in.) and Nxy = 700 kN/m (4000 lbf/in.), in which stresses are crit- 
ical for T > 0 (fig. 6(b)). The plot of final mass against the number of time 
slices is shown in figure 7 and indicates that, despite the criticality of 
stress, only three or four time slices are needed for a satisfactory answer. 


Thermal-Structural Efficiency Study 

This study is a comparative evaluation of the performance of the eight 
materials over a wide spectrum of load levels for two heating conditions. 

Except for Rene 41 and carbon/carbon, each material is protected by a layer 
of insulation. The first (low) heating condition is characterized by 
Tpeak = ”*089 K (1500° F) , and the second (high) heating condition has a value 
of Tp ea j< = 1533 K (2300° F) . In both cases T peak = sec anc3 

Tf = zOOO sec. In all calculations Ny = N x and N^ = (2/3) N x . The results 
for the low-heating condition are presented in figure 8. At high-load levels 
Gr/PI gives the lowest mass, followed by B/Al. Neither Rene 41 nor carbon/ 
carbon are competitive at high-load levels. For N x between 175 and 525 kN/m 
(1000 and 3000 lbf/in.) Rene 41 (uninsulated) and Lockalloy give the lowest 
mass, and for N x less than 175 kN/m (1000 lbf/in.) carbon/carbon gives the 
lowest mass. These results suggest possible benefits of these unprotected 
materials for lightly loaded regions, such as control surfaces. The curves for 
Lockalloy and aluminum have a corner or discontinuity that represents a change 
in the character of the final design - designs to the left of the corner are 
operating at the allowable temperature and allowable stress, whereas designs 
to the right of the corner are operating at the allowable stress but with tem- 
peratures less than the allowable value. 
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For the high-heating case in which T p ea k = 1 533 K (2300° F) (fig. 9) f 
Gr/PI again gives the lowest mass for the highest load values (above about 
1.05 MN/m (6000 lbf/in.)). For loads between 700 kN/m and 1.05 MN/m (4000 and 
6000 lbf/in.) , B/Al designs have the lowest mass. Lockalloy gives the lowest 
mass designs for N x between 262 and 700 kN/m (1500 and 4000 lbf/in.), and 
carbon/carbon gives the lowest mass for loads less than 262 kN/m (1500 lbf/in.). 
Based on these results uninsulated carbon/carbon and insulated Lockalloy again 
seem to have significant potential for application to lightly loaded regions f 
and for more highly loaded regions such as the wing, insulated Gr/PI, B/Al, 

Gr/E, and titanium seem to be worthy of consideration. These results are based 
on a simplified model and consider only strength and temperature requirements. 


Characteristics of Final Designs 

The main characteristics of selected designs in figures 8 and 9 are sum- 
marized in tables IV(a) to IV(1) . Tables IV(a) to IV(f) have results for the 
lower heating case where Tp ea ^ = 1 089 K (1 500° F) for values of N x from 
0 to -2.63 MN/m (-15 000 lbf/in.) in increments of 525 kN/m (3000 lbf/in.). 
Tables IV(g) to IV(1) have corresponding results for the higher heating case 
where Tp^^ = 1 533 K (2300° F) . The tables contain for each material the final 
mass, maximum temperature response in the final design and the time it occurs, 
and maximum stress ratio and the time it occurs. The tabulations are presented 
to provide check results for other investigators and to show some features of 
the designs. For the lowest loads, the structures operate at their allowable 
temperatures and the design problem is one of providing minimum-mass designs 
for temperature requir orients . As higher loads are considered, strength con- 
straints become important and a trade-off mechanism becomes operative, whereby 
high-temperature operation with less insulation, lower allowable stress, and 
a heavier structure is balanced against lower temperature operation with more 
insulation, higher allowable stress, and a lighter structure. For each insu- 
lated material, a load value is eventually reached where the trade-off favors 
lower temperature operation to maintain a sufficiently high allowable stress. 

As suggested by figure 3, the structural temperature in such designs is below 
the allowable temperature of the material and the designs are governed primarily 
by strength constraints. 

In sane cases, the temperature in the final design is well below the allow- 
able temperature. For example, Lockalloy under a load of N x = -2.63 MN/m 
(-15 000 lbf/in.) and ^peak = 1533 K (2300° F) attains a maximum temperature 
which is only 49 percent of its allowable temperature. This effect points out 
that contrary to conventional design practice, the minimum-mass insulated heat- 
sink structure is not necessarily one in which the structure operates at its 
allowable temperature. This observation was previously made for an ablation 
panel in reference 9. To bring out the significance of this observation, con- 
sider a set of Lockalloy designs for T = 1533 K (2300° F) at the load range 
of 0 to -2.63 MN/m (-15 000 lbf/in.), as described in table IV. A set of corre- 
sponding designs, based on the conventional practice of assuming operation at 
the allowable temperature, is generated for comparison. The latter designs are 
obtained by first sizing the structural thickness from the applied mechanical 
loads and the allowable stress at 700 K (800° F) , then fixing the structural 
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thickness, and then carrying out design calculations to size the insulation 
thickness required to maintain this temperature. Design masses from this con- 
ventional procedure and the present method are compared in figure 10. Since 
the structure has a peak temperature response at or close to the allowable tem- 
perature, there is little difference between the designs at loads below about 
525 kN/m (3000 lbf/in.). Significant differences occur at higher loads. For 
example, at 2.10 MN/m (12 000 lbf/in.) the conventional design has 40 percent 
more mass than the present design, and at 2.63 MN/m (15 000 lbf/in.) the conven- 
tional design has 47 percent more mass. These results suggest that design of 
structures for thermal applications should be based on simultaneous considera- 
tion of temperature and structural requirements. 


Effect of Heating Duration and Intensity 

This section describes a study of the effects of duration and intensity 
of heating on minimum mass. Designs are obtained corresponding to high and low 
values of intensity and duration of heating associated with winged reentry 
vehicles (ref. 18), as described by the parameters Tp ea j< and x f , respec- 
tively. Values used for Tf were 1500 sec and 2000 sec, and values for 
Tpeak were 1089 K ( 1 500° F) and 1533 K (2300° F) . 

Results of the study are presented in table V. Both intensity and duration 
of heating have significant effects on the designs. For example, increasing 
Tf for B/Al with Tp ea j c = 1533 K (2300° F) led to an increase in required mass 
of 1 8 percent. Increasing both Tp Ga ^ and Tf for the same case led to a mass 
increase of 41 percent. 


Effects of Simplifying Assumptions 

This section summarizes the effects of certain assumptions made in the 
analysis; namely, neglect of edge effects due to restraint of thermal expansion, 
neglect of time dependence of the mechanical loads, and neglect of heat loss 
from the back face of the structural layer. Details are given in appendix B. 

Restrai ned thermal expansion .- Calculations are carried out for Gr/E, 
titanium, and Lockalloy. Figure 11 shows a plot of design mass as a function 
of an edge-fixity coefficient 3 which varies between 0 and 1 . A value of 
3=1 corresponds to completely restrained expansion, and 3=0 corresponds 
to unrestrained expansion. Results indicate that the effects of restrained 
thermal expansion are negligibly small for Gr/E, are moderate for titanium (a 
22-percent increase in design mass as 3 is increased from 0 to 1 ) , and large 
for Lockalloy (an 83-percent increase as 3 is increased from 0 to 1 ) . 

Time-dependent loading .- To assess the effects of time dependence of the 
mechanical loads, design calculations are repeated for two selected cases: 

(1) Graphite/epoxy with a maximum load N x = -1.05 MN/m (-6000 lbf/in.); and 

(2) titanium with a maximum load of -1 .58 MN/m (-9000 lbf/in.) . The mechanical 
loads varied with time according to figure 12, which is representative of load 
against time during reentry (ref. 21). Designs are compared with those based 
on constant loading. Results are given in table VI, which shows that time- 
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dependent loading has a very small effect on the final design mass. Observe 
however that in the designs for time-dependent loads, the maximum stress ratio 
occurred at the times of maximum load (1600 sec) for both materials rather than 
times of peak temperature (2800 sec for Gr/E and 2000 sec for titanium) when 
loads were constant. There was no appreciable effect on the temperature 
responses for either material as a result of time-dependent loads. 

Heat transfer from back face of struc ture.- To assess the effect of 
neglecting heat transfer from the back face of the structural layer, thermal 
analyses of selected designs, which included radiative heat transfer from the 
back face, were carried out. A finite-element thermal model was formulated by 
using the SPAR thermal analyzer (ref. 22) . in the SPAR model the back face of 
the structural layer was permitted to radiate to a medium at room temperature. 
Temperature histories were calculated for the eight final designs in table IV(d) 
with the adiabatic back face and the radiating back face. Table VII contains 
a comparison of the peak temperatures for the two conditions. The largest 
change in peak temperature is 33 percent for titanium. This significant temper- 
ature difference does not necessarily translate into a large effect on design 
mass. As an example, the increase in strength of titanium from 548 K to 451 K 
(527° F to 352° F) is about 16 percent. If the structural mass of the titanium 
design is assumed to be decreased by 1 6 percent while the insulation mass 
remains unchanged, the total required mass is decreased by about 11 percent. 


CONCLUDING REMARKS 

A procedure is described for obtaining minimum-mass designs of insulated 
composite or metal structural panels subjected to transient heating and mechani- 
cal loads. The procedure uses nonlinear mathematical programming techniques, 
and analytical solutions are used for temperature histories and stresses in the 
structure. An important aspect of the present problem is the need to enforce 
constraints on transient temperatures and stresses. In the present procedure, 
the constraints are enforced at a finite number of time slices or discrete times 
over the time period of interest, and satisfactory results are obtained with a 
small number of times . 

To simplify the analysis, edge effects due to restrained thermal expansion 
are neglected, constant (time independent) mechanical loads are assumed, and 
heat loss from the back face of the structure is neglected. The errors in 
design mass introduced by these assumptions are examined. Neglecting edge 
effects can lead to sizable errors; neglecting heat loss from the back face has 
a smaller but significant effect; and time dependence of loading has a negligi- 
ble effect on the designs. 

The design procedure is used to perform structural efficiency studies for 
eight materials for two representative values of peak external heating and a 
range of mechanical load levels. Materials considered are graphite/polyimide, 
graphite/epoxy, boron/aluminum, titanium, aluminum, Rene 41 , carbon/carbon, 
and Lockalloy. All except Rene 41 and carbon/carbon have a layer of reusable 
surface insulation. For the lower heating condition graphite/polyimide gave 
the lowest mass except at the lowest load levels, whereas Rene 41 and carbon/ 
carbon gave the lowest mass. For the higher heating condition graphite/ 
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polyimide again gave the lowest mass for high and moderate loads, whereas 
Lockalloy and carbon/carbon gave the lowest mass at low-load levels • Addi- 
tional calculations indicate that intensity and duration of heating can have 
a substantial effect on design mass. 

Examinations of the characteristics of the final designs indicate that for 
sufficiently high mechanical loads, the optimum structure has a temperature 
response substantially less than the recommended allowable temperature of the 
material. This result is contrary to the usual engineering practice of basing 
a design on the higher temperature operation. Comparisons of several Lockalloy 
designs obtained by the present procedure with those based on operation at 
higher temperatures indicate that significant mass savings may be attainable by 
lower temperature operation. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 11, 1979 
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THERMAL AND STRUCTURAL ANALYSIS 


TEMPERATURE HISTORY 

The solution for the temperature history is obtained by using the analysis 
of heat transfer in a two-layered configuration given in reference 23. The pur- 
pose of this appendix is to outline the approach used in reference 23 and spe- 
cialize that general result for the present problem (fig. 1). The governing 
equations, boundary conditions, and initial conditions are 




3 2 T 

3z 2 

1 


1 


= p n C-| 


3T-) 

3t~ 


^2 


3 2 T 2 

3Z 2 " 

2 


3T 2 

P2C2 

z 3t 


3T-j 

k ' 5 ^ 


zi =0 


3T 2 

= k 2 

3z 2 


z 2 =t 2 


(Al) 


(A2) 


(A3) 


T-j ( 0, X) = T 2 (t 2 ,T) 


(A4) 


T-| ( 1 1 , T ) = T eq (X) 


(A5) 


3t 2 

3zo 


= 0 


Zo=0 


(A6) 


T - ] (z-, ,0) = T 2 (z 2 ,0) = T 0 


(A7) 


Equations (AT) and (A2) govern the temperatures in the insulation and structural 
layers, respectively. Equations (A3) and (A4) define the continuity of heat 
flux and temperature at the insulation-structure interface. Equation (A5) is 
the specification of the temperature history at the upper surface of the insula- 
tion. Equation (A6) states that no heat is lost through the back face of the 
structural layer. Equation (A7) defines the initial conditions on temperatures. 
As pointed out in reference 8, the assumption that no heat is transferred from 
the back face of the structural layer is conservative, especially when radiation 
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from the back face to internal heat sinks represents a significant heat loss. 

In such cases, however, it would probably be necessary to insulate the back face 
of the structure to protect internal equipment from excessive heating by thermal 
radiation. Further discussion of the effects of this assumption is presented 
in appendix B. Solutions to the equations are obtained in reference 23. The 
expression for the average temperature in the structural layer is of interest 
for the current problem. This solution is 


T(t) = T eq (0) + 


[T 0 - T eq (0)J ^ 
n=l 


sin Y n _ v 2 


YnP 


Yn H n 


r 

Jo 


-I 

n=l 


sin _ y 2, d _ 


Y n (P" u ) 


Yn H n 


dT eq (u) 


du 


(A8) 


du 


where y n is the nt h root of 


3£ tan Y tan - = 1 

p 


(A9) 


6 = 


k l c 2P2 t 2 
k 2 c l Pl fc l 


(A10) 


S = 


k 2tl 

k l fc 2 


(All) 


Yn 


Yn 

Yn 

H„=J- 

( K * l) 

cos y n sin — + (1 

3 

+ K) Sin Y n cos — 

P _ 


(AT 2) 


and dimensionless time is 


k 2 


P = 


2 

c 2f ) 2 t 2 


(AT 3) 
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In the present work the outer surface temperature T e g(x) 
(see fig. 2) 

is represented 

as 

T e q(x) = a e g + b e gX + c e gx2 + d e gX 3 

(X S x f ) 

(A14) 

T e q(X) = T e gf 

(x > x f ) 

(A15) 


The coefficients in equation (A14) are expressed in terms of the parameters of 
the curve in figure 2. The conditions for evaluating the coefficients are that 
the curve of equation (AT 4) passes through the three points noted on the curve 
and has a zero slope at x = Xp^j^. The results are 


T eqf T eqo + ( T peak “ T eqo) 


T f \ 


[_\ T peak' 


J eq 


T f( T peak “ T f) 


T eqo “ T peak 

c eq = ” “ 2T peak d eq 

T peak 

b e q = “^XpgajjCgg - Sx^gg^dgg 


a eq _ 


T eqo 


2 


T f 

T peak 


The coefficients b e g, c e g, and d e g are nondimensionalized by defining 

3 e q = b e g/R 

Yeq = C eq/ R ^ 

<S eq = eq/R 3 

where 


k 2 

2 

c 2P2 fc 2 


Substituting equations (AT 4) and (A15) into equation (A8) gives the following 
solutions: 
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= T eq (0) + 8 eq p + y eqP 2 + 6 eqP 3 + [t 0 - T eq (0)] ^ 


sin Y n _ v 2 


YnP 


n=l 


Yn H n 


n=l (jn ‘n 


+ 36 


eq 


?! - ^ + L(, . e -YnP) 

y 2 y 4 y 6\ 7 

» n » n 


(T ^ T f ) (Al 6) 


oo 

T = T eqf + [t 0 - T eq (0)j Y 


n=l 


e -Y^f 

*Yn H n 


-I 


sin Y n p e q 


^ Y n H n W 2 

n=l yn 

2 Y e q 


a Yn(PfP) _ e YnP 


Y, 


e -YnP _ e Yn(Pf-P)f! 


n 


- “jk e Y " (Pt ' P> 

l?n 


- r*Pf)_ 


. — e Yn (Pf“P) 

y* 


2 

+ — 
y 6 
] n 


Jn(Pf-P) _ -YnP 


(T > T f ) (AT 7) 


where 


*2 

Pf = 2 T f 

c 2P2 t 2 


In the evaluation of equations (A16) and (A17), six terms of the series have pro- 
vided sufficiently accurate temperatures. In the general case where both t-| 
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and t 2 are nonzero, the roots of equation (A9) are extracted numerically by 
a systematic search. For the special case where t-| = 0, that is, the insula- 
tion layer is omitted, the roots of equation (A9) are given by 


Y 


n 



(AT 8) 


and for this case equation (A12) reduces to 

lYir . 

- -M> n 
2/2 


H n = | n - 


(A19) 


STRESS ANALYSIS 

The stresses in the laminate are computed from elementary thermoelastic 
lamination theory (ref. 24) . The constitutive equations for a balanced symmet- 
ric laminate under inplane mechanical and thermal loading are written 


N x + N tx 


A n 

At 2 

0 



Ny + N ty 

)- 

A i 2 

A 22 

0 

< 

e y 

N X y 

V J 


0 

0 

a 66 


e xy 

L J 


(A20) 


where N x , Ny, and Njjy are constant (i.e., time independent) applied mechani- 
cal forces per unit width of the laminate and A-j -| , A-|2» A 22 / and Agg are 
laminate stiffnesses, given by 


NLAYER 



i=l 


(A21 a) 


NLAYER 



i=l 


(A21 b) 
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A 22 


NLAYER 

y Q22 fc i 

i=l 


(A21 c) 


a 66 


NLAYER 

J Qeeti 

i=1 


(A21 d) 


where is the thickness of the ith layer. The thermal forces Nj- X and 

N^-y are given by 


NLAYER 

N t x = I Q6n a x + Ql 2 a y + Ql 6 a xy ^ i ^ref) (A22) 

i=1 


NLAYER 

N ty = J [(Ql2«x + Q22«y + Q26«xy)t]i(T - T ref ) (A 23 ) 

i=l 

where T ref is the stress-free temperature and 
a x = m 2 a-| + n 2 a 2 '' 

a y = n 2 a-| + m 2 a 2 > (A24) 

a xy = 2mn(a-| - a 2 )j 

In these equations, m = cos 0, n = sin 0, 0 is the lamina ply angle, and 

Ql 1 = Qn 1 " 4 + 2 (Ql2 + 2Q66^ n2m ^ + Q22° 4 (A25a) 

Q22 = Qll n4 + 2 (Qi 2 +2Qgg)n 2 m 2 + Q 2 2 m4 (A25b) 
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Ql2 

Ql6 

Q26 

Q66 

Qll 

Q22 

Q66 

Ql2 


= (Qn + Q 22 “ ^QgeJn 2 !* 2 + Qi 2 ( n4 + m4 ) 

= (Qn - Ql2 -2Q66) 011,3 + (Ql2 ~ Q22 + 2 Q66) n3m 

“ (Qn “ Q 12 ~2Q6g)n 3 m + (Ql 2 “ Q 22 + 2Q66)ran 3 

= (Qn “ Q22 “ 2Q-| 2 - 2Q 66 )m 2 n 2 + Qgg (n 4 + m 4 ) 

E 1 > 

1 - V-] v 2 

E 2 

1 - V-j v 2 
= G 

= v iQ22 


(A25c) 

(A25d) 

(A25e) 

(A25f) 


(A26) 


where E-| , E 2 , V-j , V 2 , and G are functions of temperature. The laminate 

strains e x , ey, and exy are obtained by solving equation (A20) . The lamina 
strains are obtained as follows: 


e-j = m 2 e x + n 2 ey + mnexy 
e 2 - n 2 e x + m 2 ey - mnexy 
ei2 = 2mn(e y - e x ) + (m 2 



and the stresses in each lamina are given by 


CTi = Qn e i + Ql2 e 2 ” 
°2 ~ Ql2 e l + Q22 e 2 “ 
a 12 = Q66 e 12 


(Qll a l + Ql2 a 2) (T 
(Ql 2 a l + Q22 a 2) ( T 



(A27) 


(A28) 


The stresses in equations (A28) neglect edge effects, particularly effects of 
restrained thermal expansion due to edge supports. Estimates of the importance 
of such thermal stresses are considered in appendix B. 
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EFFECTS OF RESTRAINED THERMAL EXPANSION, TIME- DEPENDENT LOADS, AND 
BACK-FACE THERMAL BOUNDARY CONDITION 

The purpose of this appendix is to quantify effects of certain simplifi- 
cations employed in the analysis. These simplifications include neglect of edge 
effects due to restrained thermal expansion, neglect of time variations of the 
applied mechanical loads, and neglect of heat loss from the back face of the 
structural layer. 


RESTRAINED THERMAL EXPANSION 

The stress analysis used in the design procedure and described in appen- 
dix A is based on stresses and strains at a point and thus does not incorporate 
edge effects, particularly effects of thermal stresses due to restraint of ther- 
mal expansion by edge fixity. The following analysis is an attempt to account, 
in an approximate manner, for the effects on design mass of thermal stresses 
due to edge fixity. If the temperature of a panel is raised while thermal 
expansion is prevented, the panel develops forces given by 


N x = -N t x 
N y = -N ty > 
N xy = 0 


(Bl) 


when Ntx and N ty are thermal forces defined in appendix A. For the special 
case of an isotropic layer 


Eoth 

N x = N y = - (T - T ref ) 

l - V 


(B2) 


where E is Young's modulus, a is the coefficient of thermal expansion, 

V is Poisson's ratio, and h is the thickness of the layer. Equations (Bl) 
correspond to complete restraint of inplane displacement. For completely unre- 
strained expansion, the forces N x and Ny are zero. This circumstance sug- 
gests the introduction of an edge-fixity coefficient g, in which thermal 
stresses due to edge restraint are given by 


} 


(B3) 


N x = -6N tx 

Ny = -BN ty 
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where 


0 S 8 ^ 1 


Completely fixed edges correspond to 8=1, and free edges correspond to 
8=0. The analysis in appendix A is modified to incorporate equations (B3) 
by replacing N tx by (1 - 8)N tx and N ty by (1 - 8)N ty in equation (A20) . 

Evaluation of thermal-stress effects is carried out for three materials: * 

Lockalloy, titanium, and graph ite/ epoxy . These materials were chosen because 
their properties result in high, intermediate, and low values of N^x and N ky . 

(See app. A and eq. (B2) . ) The results are shown in figure 11, wherein the opti- 1 
mum mass for the three materials is plotted as a function of 8. Designs corre- 
spond to N x = -1.58 MN/m (-9000 lbf/in.), Tp eak = 1089 K (1500° F) , and 
Tf = 2000 sec. The designs for 8=0 are listed in table IV(d) . Figure 11 
indicates that thermal stress due to edge fixity may have a significant effect 
on mass. This effect is apparent for Lockalloy, wherein the mass corresponding 
to 8=1 is 83 percent higher than the mass when 8=0. The change for tita- 
nium is 22 percent and essentially zero for graphite/epoxy. The importance of 
the thermal stress term is situational dependent. These results suggest the 
need to extend the design procedure to finite panels, wherein edge effects would 
naturally be included. Another benefit from extending the current procedure to 
finite panels is to enable buckling requirements to be properly included in the 
design of insulated panels. 


TIME-DEPENDENT LOADING 

Design calculations presented and discussed in the main text are based on 
constant applied mechanical loads. In an actual application, such as reentry, 
the loads vary with time. This section describes calculations to examine 
designs of time dependence of the loads. Figure 12 depicts a representative 
plot of amplitude of force (load factor) against time during reentry (ref. 21). 
As indicated in figure 12, all three load components are assumed to have the 
same time variation. 

Calculations are carried out for two selected designs: Graphite/epoxy with 

a maximum load of N x = -1.05 MN/m (-6000 lbf/in.) and titanium with a maximum 
load of -1.58 MN/m (-9000 lbf/in.). Results given in table VI show that time- 
dependent loading has a very small effect on the final design mass. In the 
designs for time-dependent loads, the maximum stress ratio occurred at the times 
of maximum load (1600 sec) for both materials rather than at the times of peak 
temperature (2800 sec for Gr/E and 2000 sec for titanium) when loads were con- 
stant. There was no appreciable effect on the temperature responses for either 
material as a result of time-dependent loads. 


BACK-FACE THERMAL BOUNDARY CONDITION 

As stated in the main text, the analysis of the temperature history given 
in appendix A incorporates the assumption that no heat flow takes place from 
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the back face of the structural layer. This assumption follows the precedents 
set by previous analyses of similar configurations and leads to conservative 
calculated temperatures. The purpose of this section is to provide some quan- 
tification of the effects of that assumption. 

In order to carry out this quantification , a finite-element thermal model 
was formulated by using the SPAR thermal analyzer (ref. 22) . In the SPAR model, 
the back face of the structural layer was permitted to radiate to a medium at 
room temperature. Temperature histories were calculated for the eight final 
designs from table IV(d) . Table VII contains a comparison of the peak tempera- 
tures for the adiabatic and radiating boundary conditions. The largest change 
in peak temperature is 33 percent for titanium, and the smallest change is for 
Rene 41 and carbon/carbon. The effect of radiation might have been expected 
to be largest for the higher temperature materials. Certainly more heat is 
transferred from the back face of the high-temperature materials than the low- 
temperature materials. However the difference in temperatures from analyses 
with and without back-face radiation is not necessarily greater for the high- 
temperature materials. Other considerations are involved and the following 
explanation is offered. Inspection of the thermal properties and the response 
for the eight materials indicate an inverse correlation between the heat-storage 
dT 2 

term P 2 c 2 t 2 and the difference in temperatures in the two columns of 

dT 

table VII. From elementary energy conservation, a high value of heat storage 
leads to a high rate of radiation heat transfer from the back face of the struc- 
ture. This result in turn requires a high, structural wall temperature. Con- 
versely, the materials with low-heat storage (notably titanium and B/Al) tend 
to have lower wall temperatures. The temperatures in the left column of 
table VII are upper bounds to the corresponding temperatures in the right 
column. When radiation is included, the higher storage materials have struc- 
tural temperatures closer to the adiabatic values than the low-storage mater- 
ials, and hence the effect of the back wall radiation is greatest for these 
low-storage materials. 

Although the effect of the back-face boundary conditions has a fairly 
sizable effect on some peak temperatures, it does not necessarily have a large 
effect on design mass. As an example, consider the largest difference in tem- 
perature being that of titanium. The increase in allowable stress of titanium 
from 548 K to 451 K (527° F to 352° F) is about 16 percent. If the structural 
mass of the titanium design is assumed to be decreased by 16 percent and the 
insulation mass remains unchanged, the total required mass is decreased by about 
11 percent. 
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TABLE I.- THERMAL PROPERTIES USED FOR OPTIMIZATION STUDIES OF INSULATED PANELS 


Material 

k 

c 

P 


a 

W/m-°C 

Btu/in-sec-°F 

J/kg-°C 

Btu/lbm-°F 

kg/m 3 

lbm/in 3 

K 

Op 

Insulation 

0.181 

2.42 x 10 -6 

1200 

0.291 

144 

0.00521 

— 

— 

Graph i te/poly imide 

2.30 

3.08 x 1CT 5 

1297 

.310 

1 550 

.056 

533 

500 

G r aph i t e/ epoxy 

1 .64 

2.20 x 10" 5 

962 

j 

.230 

1 550 

.056 

450 

350 

Bo r on/ al urn i num 

67.8 

9.07 x 10" 4 

j 1397 

.334 

2712 

.098 

589 

600 

Aluminum 

234 

3.13 x 1(T 3 

920 

.220 

2770 

.100 

450 

350 

Titanium 

9.49 

1 .27 x 10" 4 

565 

.135 

4429 

.160 

589 

600 

Rene 41 

28.4 

1 

o 

X 

o 

00 

cn 

1260 

.302 

8249 

P 

.298 

1089 

1500 

Carbon/carbon 

6.56 

8.78 x 10 -5 

1380 

.330 

1439 

.052 

1644 

2500 

Lockalloy 

I 

167 

2.23 x 10~ 3 

I 

1 2050 

.490 

2090 

1 

.0756 

! 

700 



800 



TABLE XI.- MECHANICAL PROPERTIES OF STRUCTURAL MATERIALS USED FOR OPTIMIZATION STUDIES OF INSULATED PANELS 


Property Material 


Graphite/polyimide Graphite/epoxy Bor on/ aluminum Aluminum 


Symbol 

Unit 

RT 

533 K 
(500° F) 

RT 

450 K 
(350° F) 

RT 

589 K 
(600° F) 

RT 

450 K 
(350° F) 

Ei 

GPa 

(lbf/in 2 ) 

133 

(19.3 x 10 6 ) 

133 

(19.3 x TO 6 ) 

155 

(22.5 x 10 6 ) 

155 

(22.5 x 10 6 ) 

200 

(29 x 1 0 6 ) 

205 

(29.7 x TO 6 ) 

73.1 

(10.6 x 10 6 ) 

69.6 

(10.1 x 10 6 ) 

e 2 

GPa 

(lbf/in 2 ) 

9.10 

(1 .32 x 10 6 ) 

4.14 

(0.6 x IO 6 ) 

8.83 

(1.28 x 1 0 6 ) 

5.58 

(0.81 x IO 6 ) 

27.6 

(4 x IO 6 ) 

20.5 

(2.98 x 1 0 6 ) 

( ) 

( ) 

V 1 


0.37 

0.51 

0.30 

0.23 

0.27 

0.21 

0.33 

0.33 

G 

GPa 

(lbf/in 2 ) 

5.58 

(0.81 x io 6 ) 

4.41 

(0.64 x 10 6 ) 

5.10 

(0.74 x IO 6 ) 

0.55 

(0.08 x 1 0 6 ) 

15.5 

(2.25 x 10 6 ) 

7.03 

(1.02 x io 6 ) 

( ) 

( ) 


o c -l 

(°F _1 ) 

-0.68 x io -6 
(-0.38 x IO' 6 ) 

0.14 x 10' 6 
(0.08 x IO' 6 ) 

-0.23 x IO -6 
(-0.13 x IO -6 ) 

-0.13 x io -6 
(-0.07 x IO" 6 ) 

6.08 x 10 -6 

(3.38 x IO" 6 ) 

7.4 x IO -6 
(4.1 x 10" 6 ) 

22.8 x IO" 6 

(12.65 x IO -6 ) 

24.8 x IO -6 
(13.8 x IO -6 ) 

“2 

° C -1 

(Op-l) 

27 x IO -6 
(15 x IO -6 ) 

45 x IO" 6 
(25 x IO -6 ) 

30.4 x IO -6 
(18.9 x IO -6 ) 

78.7 x IO -6 
(43.7 x IO" 6 ) 

1.75 

(9.7 x IO -6 ) 

23.2 x IO" 6 
(12.9 x IO -6 ) 

( ) { ) 

Xrp 

GPa 

(lbf/in 2 ) 

1.09 

(1 57 400) 

1.02 

(1 47 300) 

1.11 

(1 61 000) 

1.03 

(1 50 000) 

1 .41 

(205 000) 

1 .17 

(169 000) 

0.40 
(58 000) 

0.32 
(46 400) 

x c 

MPa 

(lbf/in 2 ) 

-867 

(-1 25 800) 

-450 

(-65 200) 

-970 

1 (-1 41 000) 

-848 

(-1 23 000) 

-2530 
(-367 000) 

-1972 
(-286 000) 

( ) 

( ) 

y t 

MPa 

(lbf/in^) 

16.5 

(2390) 

6.62 

(960) 

35.8 

(5190) 

21 .0 
(3040) 

162 

(23 500) 

144 

(20 900) 

( ) 

( ) 


MPa 

(lbf/in 2 ) 

-109 

(-1 5 790) 

-87.8 
(-12 730) 

-170 

(-24 700) 

-119 

(-17 300) 

-292 

(-42 400) 

-214 

(-31 100) 

( ) 

( ) 

s 

MPa 

(lbf/in 2 ) 

93.8 
(13 600) 

53.1 

(7700) 

57.9 

(8400) 

25.9 

(3760) 

93.1 

(13 500) 

27.9 

(4050) 

( ) 

( ) 

p 

kg/m 3 

(lbm/in 3 ) 

1550 

(0.056) 

1550 

(0.056) 

1550 

(0.056) 

1550 

(0.056) 

2710 

(0.098) 

2710 

(0.098) 

2768 

(0.100) 

2768 

(0.100) 



TABLE II.- Concluded 


Property 


Titanium 


589 K 
{600° F) 


Material 


Rene 41 


1 089 K 
(1500° F) 


Carbon/carbon 


R 


1644 K 
(2500° F) 


Lockalloy 


700 K 
(800° F) 


3 21 2 117 1 3.8 20.7 1 86 

i (lbf/in 2 ) | (16.5 x 10 6 ) | (1 5 x 10 6 ) (30.8 x 10 6 ) (1 6.9 x 1 0 6 ) (2 x 10 6 ) (3 x 1 0 6 ) (27 x 10 6 ) | (16 x 10 6 ) 


(lbf/in^) 


(1 .5 x 10 6 ) | (2 x 10 6 ) 


(lbf/in^) 


.60 7.6 
(1 .1 x 10 6 ) (1 .1 x 10 6 ) 


9.5 x 10~ 6 10.4 x 10 -6 11.8 x 1 0” 6 15.2 x 1 0“ 6 3.1 x 1 0 -6 ! 3.1 x 1 0“ 6 17.6 x 1 0 -6 17.6 x 1 0“ 6 

) (5.3 x 10' 6 ) j (5.8 x 10" 6 ) (6.55 x 10 -6 ) (8.43 x 1 0 -6 ) (1.7 x 1 0” 6 ) (1 .7 x 1 0~ 6 ) . (9 .8 x 1 0” 6 ) (9.8 x 10' 6 ) 



1 

j 

3.1 x 10 -6 3.1 x 10"6 

1) ( — ) ! 

( — ) 

( — ) 

1 ( ) (1 .7 x IQ -6 ) (1 .7 x 1 0“ 6 ) 


(lbf/in^) 


(lbf/in^) 


(lbf/in^) 


(lbf/in z > ( ) 

4429 

(lbm/in J ) (0.16) 


( ) t ( ) 


(-11 500) 

(-18 300) 




4429 

8248 

8248 

i U39 

\ 

1439 

2093 

2093 

| (0.16) 1 

(0.298) 

(0.298) 

j (0.052) 1 

(0.052) 

(0.0756) 

(0.0756) 

























TABLE III.- LOWER LIMITS ON THICKNESSES USED IN CALCULATIONS 



Minimum thickness 

Material 





mm 

in. 

Insulation 

3.81 0 

0.15 

Gr aph i te/poly imi de 

.076 

*.003 

Gr aphite/epoxy 

.076 

*.003 

Boron/aluminum 

.076 

*.003 

Aluminum 

.610 

.024 

Titanium 

.610 

.024 

Rene 41 

.762 

.030 

Carbon/carbon 

.031 

*.012 

Lockalloy 

.813 

.032 


*Ply thickness 


TABLE IV.- CHARACTERISTICS OF MINIMUM-MASS INSULATED PANELS FOR VARIOUS STRUCTURAL MATERIALS 


(a) T peak = 1 089 K (1500° F) ; N x = N y = 0; N^ = 0 


Material 

Total mass 

Maximum 

temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time of 
maximum 
stress ratio, 
sec 

kg/m 2 

lbm/ft 2 

K 

°F 

G r aph i te/po ly i mi de 

11.82 

2.42 

533 

500 

2000 

0 



Gr aph i t e/epoxy 

15.87 

3.25 

450 

350 

2800 

.01 

0 

Boron/a luminum 

10.20 

2.09 

389 

600 

2000 

.10 

2000 

Aluminum 

16.06 

3.29 

450 

350 

2800 

0 

— 

Titanium 

11 .82 

2.42 

589 

600 

2000 

0 

— 

Rene** 41 

6.30 

1.29 

1078 

1480 

800 

0 

— 

Lockalloy 

7.82 

1 .60 

700 

800 

1600 

0 

— 

Carbon/carbon 

1.76 

.36 

1078 

1480 

800 

0 

— 


(b) T peak = 1089 K (1500° F) ; N x = N y = -525 kN/m (-3000 lbf/in.); 
N xy = 350 kN/m (2000 lbf/in.) 


Material 

Total mass 

Maximum 

temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time of 
maximum 
stress ratio, 
sec 

kg/m 2 

lbm/ft 2 

K 

°F 

G r aph i t e/ po ly imi de 

11.96 

2.45 

533 

500 

2000 

0.65 

2000 

Gr aph it e/epoxy 

15.87 

3.25 

450 

350 

2800 

.87 

2800 

Boron/aluminum 

11.13 

2.28 

583 

590 

1600 

1 .00 

1600 

Aluminum 

17.58 

3.60 

450 

350 

2400 

1 .00 

2400 

Titanium 

13.18 

2.70 

589 

600 

2000 

.98 

2000 

Rene* 41 

11.72 

2.40 

1078 

1480 

800 


800 

Lockalloy 

11.72 

1 2.40 

567 

560 

1600 

1.00 

1600 

Carbon/carbon 

20.41 

4.18 

1081 

1486 

800 

1 .00 

2400 
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TABLE IV.- Continued 


(c) Tpe ak = 1089 K (1500° F) j N x = N y = -1.05 MN/m (-6000 lbf/in.); 
Nxy = 700 kN/m (4000 Xbf/in.) 


Material 

Total mass 

Maximum 

temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time of 
maximum 
stress ratio, 
sec 

kg/m 2 

lbm/ft 2 

K 

°F 

G r aph i t e/poly irai de 

13,28 

2.72 

533 

500 

2000 

0.76 

2000 

Gr aph it e/epoxy 

16.60 

3.40 

450 

350 

2800 

1 .00 

2800 

Boron/aluminum 

14.84 

3.04 

522 

480 

1600 

1 .00 

1600 

Aluminum 

21 .53 

4.41 

450 

350 

2000 

1 .00 

2000 

Titanium 

16.45 

3.37 

584 

592 

2000 

1 .00 

2000 

Rene* 41 

23.39 

4.79 

1078 

1480 

800 

1 .00 

800 

Lockalloy 

16.55 

3.39 

503 

445 

1600 

1 .00 

1600 

Carbon/carbon 

41 .35 

8.47 

1076 

1477 

800 

1 .00 

2800 


(d) T peak = 1089 K (1500° F) ; N x = N y = -1.58 MN/m (-9000 lbf/in.); 
N xy =1.05 MN/m (6000 lbf/in.) 


Material 

Total mass 

Max imum 
temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time of 
maximum 
stress ratio, 
sec 

kg/ra 2 

lbm/ft 2 

K 

Op 

G r aph i te/po ly imi de 

13.96 

2.86 

530 

494 

2000 

1.00 

2000 

Graphite/epoxy 

18.02 

3.69 

450 

350 

2400 

.92 

2400 

Boron/aluminum 

16.70 

3.42 

494 

429 

2000 

1 .00 

2000 

Aluminum 

29.29 

6.00 

384 

232 

2800 

1.00 

2800 

Titanium 

20.36 

4.17 

548 

527 

2000 

1.00 

2000 

Rene* 41 

35.06 

7.18 

1078 

1480 

800 

1 .00 

800 

Lockalloy 

21 .53 

4.41 

475 

395 

1600 

1 .00 

1600 

Carbon/carbon 

62.50 

12.80 

1018 

1373 

800 

1 .00 

2800 
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TABLE IV,- Continued 


(e) Tpe ak = 1089 K (1500° F); N x = N y = -2.10 MN/m (-12 000 lbf/in.); 
*" 1 N xy = 1.40 MN/m (8000 lbf/in.) 


Material 

Total mass 

Max imum 
temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

Time of 
maximum 
stress ratio, 
sec 


kg/m^ 

lbm/ft 2 

K 

Op 

ratio 

Gr aph i te/poly irai de 

17.53 

3.59 

498 

437 

2000 

1 .00 

2000 

Gr aph i te/epoxy 

20.07 

4.11 

440 

340 

2400 

.75 

2400 

Boron/aluminum 

19.48 

3.99 

499 

438 

1600 

1 .00 

1600 

Aluminum 

33.88 

6.94 

384 

232 

2400 

1 .00 

2400 

Titanium 

24.17 

4.95 

570 

536 

1600 

1 .00 

1600 

Rene^ 41 

46.72 

9.57 

1078 

1480 

800 

1 .00 

800 

Lockalloy 

26.80 

5.49 

456 

360 

1600 

1 .00 

1600 

Carbon/carbon 

82.61 

16.92 

924 

1204 

800 

1 .00 

3600 


( f ) = 1089 K (1500° F) ; N x = N y = -2.63 MN/m (-15 000 lbf/in.); 

N^ = 1.75 MN/m (10 000 lbf/in.) 


Material 

Total mass 

Maximum 

temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time of 
maximum 
stress ratio, 
sec 

kg/m 2 

lbm/ft 2 

K 

°F 

G r aph i te/po ly imi de 

21 .73 

4.45 

417 

290 

2000 

1 .00 

2000 

Gr aph i te/epoxy 

24.90 

5.10 

450 

350 

2400 

.89 

0 

Boron/aluminum 

22.36 

4.58 

517 

370 

2000 

1 .00 

2000 

Aluminum 

38.08 

7.80 

432 

318 

2000 

1 .00 

2000 

Titanium 

28.12 

5.76 

555 

540 

1600 

1.00 

1600 

Rene* 41 

58.44 

11.97 

1078 

1480 

800 

1 .00 

800 

Lockalloy 

33.40 

6.84 

505 

449 

1600 

1 .00 

1600 

Car bon/car bon 

96.28 

19.72 

889 

1140 

1200 

1 .00 

0 



TABLE IV.- Continued 


(9> T peak “ 1533 K (2300° F) ; N x = N y = 0; Nxy = 0 





Maximum 

Time of 


Time of 


i Total mass 




Maximum 


Material 



tempe: 

rature 

maximum 
temperature , 

stress 

maximum 
stress ratio. 


kg/m 2 

lbm/ft 2 

K 

°F 

sec 

ratio 

sec 

Graph ite/polyimide 

15.77 

3.23 

533 

500 

2400 

0 

— 

G r aph i te/ epoxy 

19.53 

4.00 

450 

350 

3600 

.06 

0 

Boron/a lum inum 

13.57 

2.78 

589 

600 

2400 

.13 

2400 

Aluminum 

20.65 

4.23 

450 

350 

3600 

0 

— 

Titanium 

15.53 

3.18 

589 

600 

2400 

0 

— 

Rene' 41 

11.03 

2.26 

1089 

1500 

1200 

0 

— 

Lockalloy 

10.89 

2.23 

700 

800 

2000 

0 

— 

Car bon/ car bon 

1 .76 

.36 

1515 

2268 

800 

.96 

0 


(h) = 1533 K (2300° F) ; N x = N y = -525 kN/m (-3000 lbf/in.)? 

N xy = 350 kN/m (2000 lbf/in.) 



Total mass 

Maximum 

Time of 

Maximum 

stress 

Time of 

Material 

temperature 

maximum 

maximum 


kg/m 2 

lbm/ft 2 

K 

°F 

temperature, 

sec 

ratio 

stress ratio, 
sec 

Graph ite/polyimide 

15.77 

3.23 

533 

500 

2400 

0.26 

0 

Graph it e/ epoxy 

19.72 

4.04 

450 

350 

3600 

.93 

3600 

Boron/a luminum 

14.89 

3.05 

560 

548 

2400 

1 .00 

2400 

Aluminum 

21.43 

4.39 

450 

350 

3200 

1 .00 

3200 

Titanium 

16.40 

3.36 

589 

600 

2400 

1 .00 

2400 

Rene* 41 

15.28 

3.13 

1046 

1423 

1200 

1 .00 

1200 

Locka Hoy 

13.38 

2.74 

656 

720 

1600 

1 .00 

1600 

Carbon/carbon 

20.41 

4.18 

1520 

2277 

800 

1 .00 

0 
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TABLE IV.- Continued 


(i) Tpe ak = 1533 K (2300° F) ; N x = N y = -1.05 MN/m (-6000 lbf/in.); 
N xy = 700 kN/m (4000 lbf/in.) 


Material 

Total mass 

Maximum 

temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time of 
maximum 
stress ratio, 
sec 

kg/m 2 

lbm/ft 2 

K 

°F 

G r aph i te/poly imi de 

16.21 

3.32 

533 

500 

2400 

0.94 

2400 

Graphite/epoxy 

21 .04 

4.31 

440 

332 

4000 

1 .00 

4000 

Boron/aluminum 

16.80 

3.44 

560 

549 

2000 

.91 

2000 

Aluminum 

25.10 

5.14 

450 

350 

2800 

1.00 

2800 

Titanium 

20.94 

4.29 

589 

600 

2000 

.83 

2000 

Rene' 41 

24.17 

4.95 

713 

824 

2000 

1.00 

2000 

Lockalloy 

18.55 

u> 

00 

o 

563 

554 

2000 

1 .00 

2000 

Carbon/carbon 

41 .35 

r- 

00 

1513 

2264 

800 

1 .00 

0 


(3) Tpeak = 1 533 K (2300° F) ; N x = N y = -1.58 MN/m (-9000 lbf/in.) j 
N xy =1.05 MN/m (6000 lbf/in.) 


Material 

Total mass 

Maximum 

temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time of 
maximum 
stress ratio, 
sec 

kg/m 2 

lbm/ft 2 

K 

op 

G r aph i te/po ly imi de 

18.06 

3.70 

533 

500 

2000 

0.06 

0 

Gr aph i te/epoxy 

21 .58 

4.42 

450 

350 

3200 

.95 

3200 

Boron/aluminum 

19.24 

3.94 

516 

470 

2000 

1 .00 

2000 

Aluminum 

30.66 

6.28 

426 

307 

2800 

1 .00 

2800 

Titanium 

23.19 

4.75 

589 

600 

2000 

1 .00 

2000 

Rene" 41 

31 .64 

6.48 

783 

950 

1600 

1 .00 

1600 

Lockalloy 

23.63 

4.84 

489 

420 

2000 

1 .00 

2000 

Carbon/carbon 

62.40 

12.78 

1422 

2100 

800 

1 .00 

0 
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TABLE IV.- Concluded 


(k) T peak = 1533 K (2300° P) ; N x = N y = -2.10 MN/m (-12 000 lbf/in.); 
N xy = 1.40 MN/m (8000 lbf/in.) 


Material 

Total mass 

Maximum 

temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time of 
maximum 
stress ratio, 
sec 

kg/m 2 

lbm/ft 2 

K 

°F 

Graphite/polyimide 

19.68 

4.03 

529 

492 

2000 

1.00 

2000 

Gr aph i te/epoxy 

25.44 

5.21 

425 

305 

3600 

1.00 

0 

Boron/aluminum 

22.02 

4.51 

499 

438 

2000 

1.00 

2000 

Aluminum 

35.05 

7.18 

448 

346 

2400 

1.00 

2400 

Titanium 

27.05 

5.54 

581 

586 

2000 

1 .00 

2000 

Rene" 41 

38.96 

7.98 

709 

817 

1600 

1 .00 

1600 

Lockalloy 

28.51 

5.84 

474 

394 

2000 

1 .00 

2000 

Carbon/carbon 

82.56 

16.91 

1277 

1839 

1200 

1 .00 

0 


(1) = 1533 K (2300° P) ; N x = N y = -2.63 MN/m (-15 000 lbf/in.); 

N xy =1.75 MN/m (10 000 lbf/in.) 



Total mass 

Maximum 

Time of 

Maximum 

Time of 

Material 



temperature 

maximum 

stress 

maximum 






temperature. 

ratio 

stress ratio. 


kg/m 2 

lbm/ft 2 

K 

Op 

sec 

sec 

Graphite/polyimide 

25.00 

5.12 

454 

358 

2400 

1 .00 

2000 

Gr aph i te/epoxy 

29.34 

6.01 

433 

319 

2800 

1 .00 

0 

Boron/aluminum 

25.10 

5.14 

459 

367 

2400 

1 .00 

2400 

Aluminum 

41 .06 

8.41 

412 

282 

2400 

1.00 

2400 

Titanium 

30.71 

6.29 

532 

498 

2000 

1 .00 

2000 

Rene* 41 

49.95 

10.23 

847 

1065 

1600 

1.00 

1600 

Lockalloy 

33.79 

6.92 

463 

373 

2000 

1 .00 

2000 

Carbon/carbon 

96.28 

19.72 

1220 

1736 

1200 

1.00 

0 
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TABLE V.- OPTIMUM MASS OF INSULATED PANELS FOR VARIOUS VALUES OF 


HEATING INTENSITY. AND DURATION 

[N x = N y = -700 kN/m (-4000 lbf/in.); N xy = 525 kN/m (3000 lbf/in. fj 


Material 

T peak “ 1 

K (1500° F) 

T peak - 1 

K (2300° F) 

*■ 

Tf = 1500 sec 

Tf = 2000 sec 

x f = 1 500 sec 

T f = 2000 sec 

kg/m^ 

lbm/ft 2 

kg/m 2 

lbm/ft 2 

kg/m 2 

lbm/ft ^ 

kg/m 2 

lbm/ft 2 

Graph ite/poly imide 

11 .57 

2.37 

12.01 

2.46 

14.16 

2.90 

15.23 

3.12 

G r aph i t e/e poxy 

15.14 

3.10 

15.97 

3.27 

18.31 

3.75 

20.16 

4.13 

Bor on/a lumi num 

*13.72 

*2.81 

15.14 

3.10 

*16.36 

*3.35 

19.33 

3.96 

Aluminum 

18.16 

3.72 

19.19 

3.93 

21 .43 

4.39 

23.00 

4.71 

Titanium 

*14.26 

*2.92 

14.45 

2.96 

16.80 

3.44 

18.46 

3.78 

Rene 41 

*16.70 

*3.42 

*16.75 

*3.43 

i *- - 1 " 

— 


— 

Car bon/car bon 

*30.27 

*6.20 

*30.27 

*6.20 

*30.27 

*6.20 

*30.27 

*6.20 

Lockalloy 

*13.23 

*2.71 

*14.79 

*3.03 

*15.48 

*3.17 

*16.01 

*3.28 


♦Temperature response less than allowable value 





TABLE VI.- EFFECT OF TIME-DEPENDENT LOADS ON DESIGNS OF INSULATED PANELS 


Material 

Type of load 

Mass 

Maximum 

temperature 

Time of 
maximum 
temperature, 
sec 

Maximum 

stress 

ratio 

Time of 
maximum 
stress ratio, 
sec 

kg/m 2 

lbm/ft 2 

K 


Graphite/ 

epoxy* 

Constant 

16.60 

3.40 

450 

350 

2800 

1 .00 

2800 

Time dependent 

16.40 

3.36 

450 

350 

2800 

0.87 

1600 

Titanium** 

Constant 

20.36 

4.17 

548 

527 

2000 

1 .00 

2000 

Time dependent 

20.21 

4.14 

558 

544 

2000 

1 .00 

1600 


*N X = N y = -1.05 MN/m (-6000 lbf/in.); 
T peak = 1089 K (1500° F) . 

**N X = N y = -1.58 MN/m (-9000 lbf/in.); 
Tpeak = 1089 K (1500° F) . 


Nxy = 700 kN/m (4000 lbf/in.); T f = 2000 sec; 
Nxy = 1.05 MN/m (6000 lbf/in.); T f = 2000 sec 












TABLE VII.- EFFECT OF BACK -FACE RADIATION ON MAXIMUM STRUCTURAL TEMPERATURES 

Calculations for final designs corresponding to: 

N x = Ny = -1.58 MN/m (-9000 lbf/in.); 

N xy =1-05 MN/m (6000 lbf/in.); T f = 2000 sec; 

Tpeak = 1089 K (1500° F) 


Material 

Maximum temperature 

Adiabatic 

back face 

Radiating 

back face 

K 

°F 

K 

Op 

Graphite/polyimide 

530 

494 

453 

355 

Graph! t e/epoxy 

450 

350 

393 

248 

Boron/aluminum 

492 

429 

441 

334 

Aluminum 

384 

232 

362 

191 

Titanium 

548 

527 

451 

352 

Rene 41 

1078 

1480 

1069 

1464 

Lockalloy 

475 

395 

448 

346 

Carbon/carbon 

1 01 8 

1373 

963 

1274 
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Total 



Figure 


Optimum mass, kg/m 



Number of discrete times 

Figure 5.- Effect of number of discrete times on optimum mass of 
insulated Gr/E panel. N x = N y = -700 kN/m (-4000 lbf/in.); 
N xy = 525 kN/m (3000 lbf/in.); Tpe ak = 1089 K (1500° F) ; 

Tf = 1500 sec. 


.2 



0 1 2 3 4 5 6X 10 

Time, sec 

(a) N x = N y = -700 kN/m (-4000 lbf/in.); Nxy = 525 kN/m (3000 lbf/in.). 

Figure 6.- Temperature and stress histories for final designs of insulated 
Gr/E panel (45° ply). T nM t, = 1089 K (1 500° F) ; T f = 1 500 sec. 



Time, sec 


(b) N x = N y = -1.05 MN/m (-6000 lbf/in.); Nxy = 700 kN/m (4000 lbf/in.) 


Figure 6.- Concluded. 




0 5 10 15 20 

Number of discrete times 

Figure 7.- Effect of number of discrete times on optimum mass of 
insulated Gr/E panel. N x = N y = -1.05 MN/m (-6000 lbf/in.); 

N xy = 700 kN / m (4000 lbf/in.); = 1089 K (1500° F) ; 

Tf = 1 500 sec . 
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Figure 8.- Optimum mass of insulated panels. N y = N x ; N xy = (2/3)N x ; 

Tpeak = 1089 K (1500° F) . 
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Figure 9.- Optimum mass of insulated panels. N x = N y ; N xy 

T peak = 1533 K (2300° F) . 
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